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Abstract 

A design strategy for optimal design of composite grid-stiffened cylinders subjected to global and local buckling constraints and 
strength constraints was developed using a discrete optimizer based on a genetic algorithm. An improved smeared stiffener 
theory was used for the global analysis. Local buckling of skin segments were assessed using a Rayleigh-Ritz method that 
accounts for material anisotropy. The local buckling of stiffener segments were also assessed. Constraints on the axial membrane 
strain in the skin and stiffener segments were imposed to include strength criteria in the grid-stiffened cylinder design. Design 
variables used in this study were the axial and transverse stiffener spacings, stiffener height and thickness, skin laminate stacking 
sequence and stiffening configuration, where stiffening configuration is a design variable that indicates the combination of axial, 
transverse and diagonal stiffener in the grid-stiffened cylinder. The design optimization process was adapted to identify the best 
suited stiffening configurations and stiffener spacings for grid-stiffened composite cylinder with the length and radius of the 
cylinder, the design in-plane loads and material properties as inputs. The effect of having axial membrane strain constraints in 
the skin’ and stiffener segments in the optimization process is also studied for selected stiffening configurations. © 1998 Published 
by Elsevier Science Ltd. All rights reserved. 


1. Nomenclature 
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^sk 
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Ssk 

Su ^ 2 ? ^2 


Axial stiffener spacing 
Transverse stiffener spacing 
Stiffener height 
Skin laminate thickness 
Stiffener thickness 
Axial design load 

Design variable for stacking sequence of 
skin laminate 

Design variable for stiffening configuration 
Global buckling load factor 
Buckling load factor of skin segment 
Crippling load factor of axial, transverse 
and diagonal stiffener segment 
Strain level factor for skin segment 
Strain level factor for axial, transverse and 
diagonal stiffener segment 


2. Introduction 

An aircraft or a launch vehicle in flight is subjected 
to loads associated with all flight conditions. These 


external loads are resisted by the structure and an 
internal load distribution is established based on the 
structural layout and external loads. The internal loads, 
which depend on the location in the aerospace vehicle, 
may cause either overall buckling of the stiffened 
fuselage/fuel tank structure, buckling of the skin 
segment between stiffeners, or crippling of stiffener 
segment. A stiffened circular cylindrical shell is a 
widely used structural configuration for an aircraft 
fuselage or a launch vehicle fuel tank. An efficient and 
accurate analysis method is needed to develop a 
buckling-resistant design for the general grid-stiffened 
composite circular cylindrical shell and to identify the 
most effective stiffening configurations and stiffener 
spacings for shells subjected to combined in-plane 
loading conditions. The identification of structurally 
efficient stiffening configurations and stiffener spacings 
also requires integration of optimization techniques 
with accurate structural analysis methods. 

The optimum design of stiffened panels or shells 
that satisfy buckling constraints either with or without 
strain constraints has been of considerable interest to 
researchers since aircraft and launch vehicle structures 
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consist mostly of stiffened panels or shells. Researchers 
have used gradient-based optimizer, discrete optimizer, 
and the simplex method. A gradient-based optimizer 
were used in Refs [1-8], a discrete optimizer was used 
in Refs [9-13], and the simplex method was used in 
Refs 1 14], [15]. 

A literature survey of the optimal design of stiffened 
panels and shells is presented in Ref. [13]. The survey 
deals with the type of optimizers that have been used, 
and the type of buckling analysis, and stiffened panel 
or shell optimized. Type of stiffened panel or shell 
refers to whether the panel or shell is stiffened axially, 
orthogonally or in multiple directions (e.g. grid- 
stiffened panel or shell). According to the review, 
optimization of panels or shells stiffened in multiple 
directions were considered in Refs [5], [6], [8], [16], 
Grid-stiffened cylinders were considered in Refs [6], 
[8], [16], and geodesically stiffened flat panels were 
considered in Ret. [5], More recently grid-stiffened 
curved panels were considered in Ref. [13], 

The survey also shows that stiffener spacing is 
treated as a design variable in Refs [I], [2], [7], [8], 
[14], [15], A gradient-based optimizer was used in Refs 
[1], [2], [7], [8], whereas the simplex method [17] was 
used in Refs [14], [15]. The stiffener panel spacing is a 
discrete variable since it can be a certain multiple of 
the length or width of the panel or shell. Both the 
gradient-based optimizer and the simplex method 
assume continuous variables and in addition the 
gradient-based optimizers require derivative informa- 
tion. Therefore as discussed in Jaunky et al. [13] 
rigorous optimization with respect to stiffener spacing 
as a design variable is best achieved using a discrete 
optimizer such as the genetic algorithm [II], 

Reddy et al. [16] studied the buckling response of 
isogrid and orthogrid cylinders with various stiffener 
spacing using global and local analyses and concluded 
that, in general, isogrid cylinders are more efficient 
than orthogrid cylinders. It is also shown [7] that 
optimizing axially stiffened panels for different fixed 
stiffener spacings using a gradient-based optimizer can 
be a good strategy. For panels or shells stiffened in 
multiple directions, this approach is tedious since one 
has to consider many starting points in the optimal 
design process. The stiffener spacings and stiffening 
configurations can be treated as design variables in an 
optimization process provided a discrete optimizer is 
used. A discrete optimizer such as the genetic 
algorithm has been used to optimize unstiffened or 
axially stiffened composite panels [9-13], and for other 
engineering problems [18], Haftka et al. [9-12] showed 
the efficiency of the genetic algorithm in dealing with 
discrete design variables for optimizing composite 
panels. They treated the skin and/or stiffener laminate 
as discrete variables since in practical applications, the 
ply orientations are limited to 0, 90 or +45° and the 


lamimtc thickness can only be integer multiples of 
comm mcially available ply thickness. Jaunky et al. [13] 
optimized flat and curved grid-stiffened panels with 
stiffen ;r height and thickness, skin laminate, stiffener 
spacings and stiffening configuration as design 
variables with a global buckling constraint only and 
showed the efficiency of the genetic algorithm in 
dealing with all these discrete design variables that 
affect he buckling response of grid-stiffened panels. 

The literature survey [13] also comments on the 
bucklii g analyses that have been used. Finite element 
analysis [6], discrete analysis [5], PASCO [19] type 
analysis [3], and global and local analyses. The merits 
of these types of analyses are discussed in Jaunky et al. 
[13], Most researchers have used global and local 
analyses, for example Refs [1], [2], [7], [8], [13-16], 
Global and local analyses can be made computationally 
very efficient. Researchers have neglected the aniso- 
tropic properties of skin segments in the local buckling 
analysis and also the curvature of the skin segments. 
Restrictions on the geometry of the skin segment and 
loading were imposed in order to obtain a closed-form 
solution for the buckling response. In some cases the 
eccentricities of the stiffener were even neglected. The 
main problem of assessing the local buckling of the 
skin segment was the non-availability of computation- 
ally efficient and accurate buckling analyses for curved 
skin segments with anisotropic material properties, and 
differer t planform geometry (e.g. general triangular 
and arf itrary quadrilateral geometry). In Jaunky et al. 
[13] accurate buckling analyses [21—23] were used for 
global and local analyses, and the occurrence of 
non-sirr ultaneous global and local failures were used as 
this condition leads to designs that are less sensitive to 
imperfections as discussed in Simitses [15], 

The present paper presents the analysis and weight 
optimization strategy for grid-stiffened composite 
circular cylindrical shells subjected to axial load and a 
global buckling design constraint as well as strength 
constrains using the genetic algorithm. Design 
variable; are the stiffener height and thickness, skin 
laminate , stiffener spacings and stiffening configura- 
tion. It is intended to be a design tool that can be used 
as a preliminary design stage for grid-stiffened cylin- 
ders for aircraft fuselage or launch vehicle application. 

In the authors’ opinion such a tool for cylinders 
stiffenec in multiple directions is lacking. The results 
are presented and discussed for cylinder designs with 
and without strength or strain constraints which 
explains the damage tolerant characteristics of grid- 
stiffened structures [20], 

2. 1. Sheil buckling analysis 

The 1 uckling analysis of grid-stiffened composite 
shells su >jected to combined loads requires several key 
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steps. In the present study, acceptable designs are 
those which buckle globally and do not exhibit any 
local skin buckling or stiffener crippling. The first step 
in the design process is to assess the global buckling 
response of a grid-stiffened shell. Once this global 
buckling response is determined, the second step is to 
determine the local skin buckling response for the 
quadrilateral and/or triangular skin segments between 
the stiffeners. The third step is to determine whether 
stiffener buckling or stiffener crippling has occurred at 
this global buckling load level. 

The global buckling analysis is based on a Rayleigh- 
Ritz method using a first-order shear-deformation 
theory and the improved smeared-stiffener modeling 
approach discussed in Ref. [21] that accounts for skin- 
stiffener interactions. It is shown [21] that buckling 
loads for stiffened panels with different stiffening 
configurations obtained using the improved smeared 
stiffener theory are in good agreement with results 
obtained using detailed finite element analysis. The 
cylinder is assumed to be simply supported and hence 
the Rayleigh-Ritz method for the global analysis 
assumes the following mode shape for the transverse 
deflection w: 

N ( n t y ny \ . m,nx 

w=Zl /W sin — 4- B mtli cos — Ism — — ( 1 ) 


where L and R are the length and radius of the cylin- 
drical shell, respectively, and N is the number of terms 
in the Fourier series. The coordinate system for the 
cylinder is shown in Fig. 1. 

The buckling analysis of local skin segments is also 
based on a Rayleigh-Ritz analysis using a first-order, 
shear-deformation theory and accounts for material 
anistropy. Boundary restraints on the skin segments are 
provided by the stiffeners and hence, the analysis must 
be capable of accommodating a variety of boundary 
conditions and a variety of skin-segment shapes [22], 
[23]. In most cases, the skin segments for grid-stiffened 
panels will have either a general parallelogram- or a 
general triangular-shaped planform. The skin segments 
are assumed to have simply supported boundary condi- 
tions in the present study. It is shown [22], [23] that the 
buckling analyses presented for arbitrary anisotropic 
quadrilateral plates and general triangular anisotropic 
plates can accommodate different boundary conditions, 
and provide buckling loads that are in very good agree- 
ment with finite element analysis and existing solutions. 
The Rayleigh-Ritz buckling analysis method for global 
and local analysis makes use of the Sanders-Koiter 
shell theory [24], [25]. This shell theory provides 
buckling loads that are in good agreement with finite 
element analysis compared to other shell theories as 
shown in Refs [26], [27]. 


STIFFENED 

CYLINDER 



Fig. 1. Unit cell of a grid-stiffened cylinder showing design variables. 
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In addition to analyzing the local skin segment for 
buckling, the local stiffener segments must be analyzed 
to determine whether stiffener crippling will occur. 
Reference [16] provides a method for determining the 
crippling load of a stiffener segment. Accordingly, the 
stiffener segment is assumed to be clamped at the 
nodes or intersection points of the stiffeners while the 
edge of the stiffener along the stiffener-skin attach- 
ment line is assumed to be simply supported. 

The global buckling load is assumed to be a scalar 
multiple of the design load and has the form 

*, = Ag Ni (2) 

where A, is the applied prebuckling axial load and 
represents the design load. Once the global buckling 
load factor (A G ) has been determined using the 
improved smeared stiffener theory, the loads acting on 
the stiffener and skin segments have to be determined 
by distributing the loads based on the extensional stiff- 
ness of the skin and the stiffener. The procedure for 
distributing the applied loads for a general grid- 
stiffened panel is discussed in Refs [5], [26] and is the 
same for a general grid-stiffened circular cylindrical 
shell. The loads acting on the skin and stiffener 
segments are computed based on a global load of 
N x = ?. G N l and these loads are used to determine the 
local buckling load factor of the skin (2 sk ), local 
crippling factors of axial stiffener segment (/,), trans- 
verse stiffener segment (A 2 ) and diagonal stiffener 
segment, (/ 3 ). 

critical load of local segment 

/,= : — — (3) 

load in local segment due to 2 G xN } 

where /= 1, 2, 3 and sk. These local buckling and 
crippling load factors describe the buckling character- 
istics of the stiffened cylinder and is as follows 

• For A sk , a 2 , x 3 >1.0, the cylinder buckles globally 
at an axial load of AqN\, i.e. A cr = a c . 

• If one of x sk , x h 2 2 , 2 3 < 1 .0, then the stiffened 
cylinder buckles locally. If x sk <1.0, then skin 
buckling occurs and if A t <1.0 then crippling of the 
axial stiffener occurs. For this case, x cr = i,x A G 
where a, <1.0, and i can be any one of sk, 1,2 or 3. 

• If more than one of x sk , 2,, x 2 and x 3 are <1.0, then 
local buckling of the stiffened cylinder occurs and 
/ tr = A,x A c , where is the minimum of any of x sk , a u 
a 2 or / 3 with values <1.0. 


3. Strain analysis 

The critical buckling load of the stiffened cylinder is 
/ cr iV t where x cr takes on values as discussed previously 
and based on this load value the loads acting on the 
skin and stiffener segments are obtained [5], [26]. For 


an axial load in the skin segment of N^ k , and the loads 
in axial, transverse and diagonal stiffener segments of 
N x u N. 2 and N x3 , respectively, the axial membrane strain 
the skm and stiffener segments are 

4sk = & M k) A^ rsk + fl ( i2 k W >sk + « | S 6 k W xysk (4) 

2 

6t| 3 iW jr3 

where 4*, 4t> 4 2 and 4* are the axial membrane strains 
in the skin segment, axial, transverse and diagonal stiff- 
ener segments respectively. The quantities a[* k) , a\\\ a$ 
and a <3) are axial flexibilities of the skin, axial, trans- 
verse and diagonal stiffeners respectively. 

The strain level factors for the skin, axial, transverse 
and diagonal stiffener segment are 

5 S k=(4skL/4sk (5) 

Sl = (4st)a/4l 

s 2 = (4A//4 2 

S 3 “ (4st)a/43 

where (e Xfik ) al and are the allowable axial 

memb ane strains in the skin and stiffeners, respec- 
tively. The values for (4 S k)u/ and (e^ t ) a/ are taken from 
D. R. Ambur (pers. commun.). 

These global and local buckling analysis methods 
and the strain analysis have been integrated into a 
computer code to provide a computationally efficient 
tool for predicting the buckling load and the strain 
level factors for grid-stiffened composite circular cylin- 
drical ihells subjected to axial compression. 


4. Shell design procedure 

The design variables for a grid-stiffened composite 
shell are the axial and transverse stiffener spacings (a, 
b ), the stiffening configuration (ICON), which is the 
combination of axial, transverse and diagonal stiff- 
eners, the skin laminate (LAMI), and the height (h) 
and thickness (t s ) of the stiffener (see Fig. 1). Except 
for the height of the stiffener, these design variables 
take 01 discrete values. Hence gradient-based optimi- 
zation methods are not suitable for optimal design 
grid-stiffened shells. Furthermore, owing to manufac- 
turing constraints a ‘family’ of good designs is needed 
rather than a single point design. 

The genetic algorithm is a method for ‘evolving’ a 
given iesign problem to a family of near-optimum 
design? [10], [11], [13], [18]. Stochastic processes are 
used O generate an initial population of individual 
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designs and the process then applies principles of 
natural selection and survival of the fittest to find 
improved designs. Furthermore, since the discrete 
design procedure works with a population of designs it 
can explore a large design space and climb different 
hills. This is a major advantage as the converged 
solution may contain many optima of comparable 
performance. The cost of having a large number of 
function evaluations is offset by the fact that a large 
number of optimum solutions are now available. The 
population or family of good designs produced by using 
the genetic algorithm may include the global optimum 
design, rather than a single design. Hence, it is an 
appropriate tool for designing general grid-stiffened 
composite shells. 


5. Design problem definition 

The present design problem is to minimize the 
weight of a grid-stiffened composite circular cylindrical 
shell given the design loading condition, the length and 
radius of the cylinder and the material properties for 
the skin and stiffeners. The design variables include 
stiffener spacings (a, b), the stacking sequence of the 
skin, stiffener layout, stiffener thickness (/„), and stiff- 
ener height ( h ) as shown in Fig. 1. All stiffeners are 
assumed to be of the same height and thickness for 
manufacturing and assembly purposes. The design 
sought here is a cylinder of minimum weight in a 
certain design space with buckles globally at the design 
loads while the axial membrane strain in the skin and 
the stiffener segments do not exceed the allowable 
axial membrane strain kb , 1 and respectively. 

This design problem cam be defined by setting up the 
optimization procedures in the following way. First, the 
global buckling load is assumed to be a scalar multiple 
of design loads and has the form 

n x =a g n„ < 6) 

where N, is the applied in-plane prebuckling load. This 
value represents the design loads for the grid-stiffened 
cylinder. Second, the design constraints imposed on 
panel include: 

• The critical buckling load should be greater than or 
equal to the design loads, that is, A c > 1 . 

• Skin segments should not buckle at the critical 

buckling load, that is, /*>!. 

• Stiffener segments should not cripple at the critical 

buckling load, that is, Ai, A 2 , A 3 ^ 1 • 

• The axial membrane strain in the skin segment 

should be less than or equal to that is, 5^^ 1 . 

• The axial membrane strain in the stiffener segment 
should be less than or equal to (&)«/, that is, 5,, S 2 , 
S 3 > 1. 


The general form of each constraint equation is written 
as 


' (l «'-" sa0 ,.l N, 

(l/Sj— 1 ) < 0.0 


(7) 


Finally, the ‘fitness’ expression based on exterior 
penalty function approach is 


Fitness = 


Q 

F(\,r ,j 


= Max 


Q 

W(X) + r, t [|g/X)| + g/X )] 2 


( 8 ) 


• where X = design variable vector 

• F(X,r,) = modified objective function 

• IF(X) = weight of panel per unit area 

• r, I [|ft(X)| +g,(X)] 2 = penalty function 

) 

. <2 = normalizing constant 

• /V c — number of design constraints 

• r, = penalty parameter 

• i = generation or iteration cycle in the optimization 
procedure. 


6. Design process based on genetic algorithm 

Implementation of the genetic algorithm [11] is 
shown schematically in Fig. 2. The design process 



Fig. 2. Flow chart for the optimization procedure using the genetic 
algorithm. 
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begins with a random selection of a specified number 
of designs which comprise the initial population (i.e. 
first generation) for the genetic algorithm. Material 
properties, radius and length of the cylinder, boundary 
conditions of the skin segment and design loadings are 
input to the analysis processor routine. The buckling 
analysis is performed which provides the critical eigen- 
values for the global buckling response of the grid- 
stiffened cylinder, the local buckling response of the 
skin and stiffener segments, and the strain level factors 
of the skin and stiffener segments. The weight of the 
grid-stiffened cylinder is also computed. This 
procedure is repeated for each design configuration in 
the population. The ‘fitness' processor then evaluates 
the fitness of each design using eqn (8) and assigns a 
rank based on the fitness expression or objective 
function. The current population of design configura- 
tions is then processed by the genetic operators (cross- 
over, mutation and permutation) to create a new 
population of design configurations for the next gener- 
ations which combines the most desirable character- 
istics of previous generations. Designs from previous 
generations may be replaced by new ones (i.e. 
children) except for the ‘most fit' designs (i.e. parents) 
which are always included in the next generation. The 
process is repeated until design convergence is 
obtained, which is defined herein by specifying a 
maximum number of generations (NSTOP) that may 
occur without improvement in the best design. The 
design procedure will now be demonstrated on grid- 
stiffened cylinders with and without the axial 
membrane strain constraint in the skin and stiffener 
segments to study the influence of constraints on the 
axial membrane strain on the optimal designs. 

6. 7. Numerical results for grid-stiffened cylinders 

Results are presented for composite grid-stiffened 
cylinders subjected to axial compression in order to 
demonstrate the changes in optimal designs due to 
constraints on the axial membrane strain. These 
changes will depend on the loading, stiffening 
configuration ( ICON ), stiffener spacings ( a , b) and the 
skin laminate (LAMI). The cylinders studied in 
examples 1 and 2 are 7391 mm long, have a radius of 
2426 mm and have blade stiffeners made of unidirec- 
tional materials [16]. The nominal ply mechanical 
properties are: E u = 63218 N mm 2 ; £\ 2 = 5604 

Table 2 

Design space for a. b, h and /, for Examples 1 and 2 

211.2 205.3 199.7 

508.0 491.6 476.3 

10.46 10.80 11.11 

1-37 1.52 1.68 


a, mm 217.4 

b, mm 525.5 

h, mm 10.16 

/„ mm 1.22 


Table 1 

Design >pace for design variable ICON and LAMI 

Integer 

LAMI 

ICON 

value 

1 

1 ± 45/0K 

Axial stiffeners 

2 

[±45/90]* 

Axial stiffeners 

3 

[±45/0/90]* 

Axial and transverse stiffeners 

4 

| ± 45/0 2 ] 2s 

Diagonal stiffeners 

5 

[±45/90,]* 

Axial and diagonal stiffeners 

6 

[ F45/02/90J* 

Transverse and diagonal stiffeners 

7 

[±45/0/90 2 ]* 

Axial, transverse and diagonal stiffeners 

8 

[±45/02/90;]* 

No stiffeners 

N mm 

, G\2 = G|j 

= G 2 3 = 2973 N mm' 2 and 

v’i 2 0.300. The mass density of the material p is 

1.578 F-06 kg mm 5 . 

The probabilities used for cross- 


over, mutation and permutation are 1.0, 0.10 and 0.95, 
respectively, and the number of discrete values that 
each design variable can accommodate is eight, which 
is a limitation of the FORTRAN code for the genetic 
algorithm being used. The design variables ICON and 
LAMI are described in Table 1 

6.2. Example 1 

The irst example is a cylinder subjected to an axial 
load of N, = 175.2 N mm 1 which represents a cylinder 
case stidied in Ref. [16], The design variables are the 
axial ard transverse stiffener spacings, the height and 
thickness of the stiffeners, the stiffening configuration 
(ICON) and the skin laminate (LAMI). The thickness 
of each ply of the skin laminate for this example is 
0.1524 mm. The design space for a , b, h and t s is shown 
m Table 2. The height and thickness of the stiffener 
(h/t) is chosen such that the aspect ratio of the stiffener 
is between 4 and 10, and the stiffener thickness is an 
integer multiple of 0.1524 mm due to manufacturing 
constraints. The minimum values of the stiffener 
spacings are also due to manufacturing constraints. A 
population size of 20 is used, NSTOP = 15, and the 
penalty parameter r, = 100000 at any iteration. The 
aliowabl : strains (i:,^) al and are 2428 microstrain 
and 1092 microstrain, respectively (D. R. Ambur, pers. 
commun). 

The lesults from design optimization with and 
without strain constants are shown in Table 3. The 
minimum weight for this cylinder is 481.53 kg 
compare! to the weights of 479.90kg for an isogrid 


!94.5 189.5 

161.8 448.2 

1143 11.75 

1.83 1.98 


184.8 180.3 

435.4 423.3 

12.07 12.38 

2.13 2.29 
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Best designs obtained by genetic algorithm for grid-stiffened cylin- 
ders (N x = 175.2 N mm l ) 


Best designs obtained by genetic algorithm for grid-stiffened cylinder 
with no constraints on strain and ICON = 5, (N % = 315.3 N mm ) 


Design variables (mm) 

Weight (kg) 

Aj 

Si 

Design variables 
(mm) 

Weight 

(kg) 

/„j 

a = 189.5, b = 525.5 
h = 10.16, 6 = 1.22 
LAMI = 4, ICON = 6 

481.53 

A g = 1.008 
A* = 1.009 
Ai = 55.1 

5* = 2.68 
S 3 = 105.1 

a = 217.4, b = 423.3 
h - 11.43, 6 = 1.83 
LAMI = 4 

634.1 

= 1.0088 
2* = 1.078 
2, = 1.11, 2, = 24.6 

a = 189.5, b = 491.6 
h = 10.16, r s = 1.22 
LAMI = 4, ICON = 6 

481.72 

Ao= 1.025 
A* =1.011 
A, = 43.4 

S sk = 2.63 
$3 = 81.7 

a = 194.5, b = 435.4 
h = 10.16,6 = 1.98 
LA IM/ = 3 

634.8 

2., = 1.0447 

2.. = 1.051 

2, = 1.04, 2., = 35.3 

a = 189.5, b = 476.3 
h = 10.16, 6 = 1.22 
LAMI = 4, ICON = 6 

481.5 

A c; = 1.034 
A* =1.012 
A 3 = 38.8 

S sk = 2.61 
5, = 72.6 

a = 194.5, h = 423.3 
h = 10.16, 6 = 1.98 
LAMI = 3 

635.5 

2,,= 1.0541 
2* =1.077 
2, = 1.04, 2, = 32.1 


cylinder and 477.63 kg for an orthogrid cylinder [16]. 
Only three optimal designs are identified here and they 
all buckle globally at their respective A G values since A sk 
and X) values are greater than 1. For this case, the 
optimal designs have ICON = 6, that is the stiffening 
configuration has transverse and diagonal stiffeners 
only, and LAMI = 4, which corresponds to a skin 
laminate which a stacking sequence of [+45/0 2 ] 2s , 
h = 10.16 mm and / s = 1.22 mm. In this case, the 
optimal designs without strain constraints are the same 
as the optimal design with strain constraints. The 
buckling-resistant design for this load case results in 
in-plane stiffness values that provide large strain level 
factors. This analysis suggests that the cylinder design 
satisfies both buckling and strength requirement. 
Convergence is obtained after 49 generations for the 
case without strain constraints and after 77 generations 
for the case with strain constraints. 

6.3. Example 2 

The second example is for a grid-stiffened cylinder 
as in Example 1 but subjected to an axial load of 
Ai, = 315.3 N mm '. For this case, the stiffening 
configuration with axial and diagonal stiffeners ( ICON) 
only has been selected. The design variables are the 
axial and transverse stiffener spacings, the height and 
thickness of the stiffeners, and the skin laminate with 
each ply being 0.2032 mm thick. The design space for 
the skin laminate is the same as shown in Table 1, and 
the design space for a, b, h and / s is shown in Table 2. 
The material properties, the genetic parameters and 
the strain allowables for Example 2 are the same as in 
Example 1. The optimal designs without strain 
constraints are shown in Table 4, whereas the optimal 
designs with strain constraints arc shown in Table 5. 
Only three optimal designs are shown here for each 
case. 

All the designs in Tables 4 and 5 buckle globally 
since /l sk , Ai and A 3 are all greater than one. The 


optimal designs obtained with strain constraints are 
^4-6 kg heavier than the optimal designs obtained 
without strain constraints for approximately the same 
global buckling load. The designs in Table 5 have 
larger axial stiffener spacing (a) than the second and 
third designs in Table 4. For the designs in Table 5, the 
strain level factors of the axial stiffeners (5,), are very 
close to unity whereas the strain factors for the 
diagonal stiffeners (5,) are much larger. The increase 
in structural weight of the optimal design for this case 
over the designs in Example 1 is due to the increase in 
load and the choice of stiffening configuration. Even 
for this load condition, imposing strain constraints docs 
not result in a significant increase in weight. Converg- 
ence is obtained after 18 generations for the case 
without strain constraints and after 41 generations for 
the case with strain constraints. 

6.4. Example 3 

The third example is a grid-stiffened cylinder 
subjected to an axial load of N\ = 350.4 N mm \ This 
represents a generic, wide-body transport aircraft 


Best designs obtained by genetic algorithm for grid-stiffened cylinder 
with constraints on strain and ICON = 5 (N x = 315.3 N mm ) 


Design variables 
(mm) 

Weight 

(kg) 

Ai 

S { 

a = 21 1.2, b = 423.3 
h= 11.11,6= 198 
LAMI ' = 4 

638.0 

Ati= 1.0113 

A* = U01 

A| = 1.36, A 3 = 32.9 

S sk = 2.07 

Si = 1-08 

S 3 = 27.2 

a = 205.3, h = 423.3 
h= 11.11,6= 198 
LAMI = 4 

639.0 

A Ci = 1-0100 

A* =1.126 

Ai = 1.38, Ai = 36.2 

$* = 2.08 
$,= 1.08 
$, = 29.8 

a = 21 1.2. b = 423.3 
h = 1 1.75, 6= 1.98 
L4MI = 4 

641.5 

Ao = 1.0600 

A* =1.053 

A, = 1.19, A 3 = 28.9 

$*=1-99 
$,= 1.03 
$ 3 = 26.0 
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fuselage section between two frames. The radius of the 
cylinder is 2426 mm and the length is 559 mm. The 
design variables are the transverse stiffener spacing, 
the height and thickness of the stiffeners and the skin 
laminate with each ply being 0.2286 mm thick. The 
optimization is performed with the stiffening configura- 
tion selected to ICON = 1 and ICON = 5. For the case 
when ICON = 5, the axial stiffener spacing is 
159.7 mm. The design space for the skin laminate is the 
same as shown in Table 1 and the design space for b. It 
and / s is shown in Table 6. The material properties and 
the genetic parameters for Example 3 is the same as in 
Example I. Three optimal designs obtained with strain 
constraints are shown in Table 7 for the axially 
stiffened cylinders and in Table 8 for the cylinders 
stiffened with axial and diagonal stiffeners. For axially 
stiffened cylinders convergence is obtained at 38 gener- 


ations. while for the cylinders with axial and diagonal 
stiffen jrs, convergence is obtained at 20 generations. 

All the optimal designs in Tables 7 and 8 buckle 
globally and do not violate any strain constraint. 
Howe\er, the optimal designs with ICON = 5 are 
slightly heavier than the optimal designs with 
ICON = 1. The buckling behavior of the first design in 
Table 7 and Table 8 subjected to an axial load of 
NX = 350.4 N mm 1 and combined transverse compres- 
sion (N y ), and shear (N^) is sought to determine their 
capability to support additional in-plane loads. A load 
combination for N x = 350.4 N mm “ N y and N v is 
sought in such a way that the cylinder can still support 
an axial load of N x = 350.4 N mm' 1 , i.e. 4 r >l and 
/c r a:l. The results are shown in Table 9 for the axially 
stiffened cylinder and in Table 10 for the grid-stittened 
cylinder with ICON = 5. It is seen that the grid- 


Table 6 


Design space for b , /i, r s for Example 3 


Axially stiffened cylinder, ICON = 1 


b (mm) 
h (mm) 
t s (mm) 

217.4 

16.19 

2.59 

211.2 

16.51 

2.74 

205.3 

16.83 

2.89 

199.7 

17.14 

3.05 

194.5 

17.46 

3.20 

189.5 

17.78 

3,35 

184.9 

18.10 

3.50 

180.3 

18.42 

3.66 




Grid-stiffened cylinder, ICOh 

= 5 




b (mm) 
h (mm) 
(mm) 

324.3 

13.36 

1.68 

319.0 

13.65 

1.83 

311.0 

13.97 

1.98 

304.8 

14.29 

2.13 

298.8 

14.67 

2.29 

293.1 

14.92 

2.44 

287.6 

15.24 

2.59 

282.2 

15.56 

2.74 


Table 7 

Best designs obtained by genetic algorithm for axially cylinder with constraints on stra n, (N, = 350.4 N mm-') 


Design variables (mm) 

h = 217.4, h = 16.51 
I, = 2.59, LAMl = 4 
b = 211.2, h = 16.19 
/, = 2.59, LAMl = 4 
b = 21 1.2, h= 16.51 
/, = 2.59. 1A Ml = 4 


Weight (kg) 
54.4 
54.4 
54.6 



Xo = 1.0177 

4k = 2.803, 2, = 1.025 

2g = 1.0079 

4* = 3.004, 2, = 1.073 

4j = 1.0285 

4k = 2.952, 2, = 1.018 


S lk = 2.310 
S,= 1.200 
S A = 2.335 
S, = 1.213 
S«k = 2.295 
5, = 1.192 


I able 8 

Best designs obtained by genetic algorithm for grid-stiffened cylinder with constraints c n strain, ICON = 5 (N, = 350.4 N mm) 


Design variables (mm) 

Weight (kg) 

b = 319.0, h = 13.36 
/,= 1.98, LAMl = 4 

55.8 

b = 311.0, h= 13.36 
6= 1.98, LAMl = 4 

55.9 

b = 304.8, h = 13.36 
/,= 1.98. LAMl = 4 

56.0 




- o - 1.0368, 4k = 2.341 S jk = 2.08 

- 1 ~ 1 .068, 2j = 24.6 Si = 1 98 

S, = 27.0 

/ u - 1 .0490, 4, = 2.387 £ sk = 2.064 

/ 1 = 1 .057, 2, = 22.79 5, _ j 072 

S, = 24.67 

/ - 1 .0595, 4k = 2.424 5 sk = 2.048 

/,= 1.049, 2, = 21.44 5, = 1.064 

S, = 22.96 
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Table 9 

Buckling load factors and strain level factors for axially stiffened cylinder 


Loading (N mm ') N x = 350.4 N x — 350.4 N x — 350.4 

N y = 3.5 N y = 2.6 N v = 0.9 

= 0 N v = 9.6 N„ = 28.9 


2 , 

Si 


/lo = 1.0085 
2* k = 2.710 
2! = 1.034 
^ = 2.349 
5, = 1.211 


2a = 1.0108 

2, k = 2.725 
2, = 1.032 
S sk = 2.339 
5, = 1.208 


2g= 1.108 

2 sk = 2.757 
2, = 1.032 
S sk = 2.330 
5, = 1.208 


= 350.4 
Ny = 0 
^> = 38.5 


2o= 1.0087 
2^ = 2.227 
2, = 1.034 
S sk = 2.330 

Si = 1.211 


stiffened cylinder can support up to 70.1 N mm 1 of 
transverse compression and 70.1 N mm -1 in shear 
whereas the axially stiffened cylinder can only support 
additional loads of 3.5 N mm -1 in transverse compres- 
sion and 38.5 N mm -1 in shear without exceeding the 
strain allowables. It is assumed here that the strains ey 
and jxy, due to N y and do not produce strains large 
enough to create failure since N y and are small 
compared to N x . For these loading combinations, both 
cylinders buckle globally and the constraint on the axial 
membrane strains in the skin and stiffener are not 
violated. Hence the grid-stiffened cylinder can sustain 
larger magnitudes of additional transverse compression 
and shear loads than the axially stiffened cylinder. 


without strain constraints for simply supported cylin- 
ders that buckle globally. It is also shown that a grid- 
stiffened cylinder optimized for axial load has an 
additional load capacity for combined loading larger 
than for a conventional axially stiffened cylinder. This 
additional load capability permits load redistribution in 
the event of structural damage and suggests that a grid- 
stiffened cylinder is more damage tolerant than an 
axially stiffened cylinder. 
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7. Concluding remarks 

A minimum-weight design optimization tool with 
global buckling constraint and constraint on the axial 
membrane strain in the skin and stiffener segments has 
been developed for grid-stiffened cylinders using global 
and local buckling analyses and a genetic algorithm. 
Design variables used are axial and transverse stiffener 
spacing, stiffener height and thickness, skin-laminate 
stacking sequence and stiffening configuration. The 
present design optimization procedure has been 
validated by comparing with existing results. Results 
for grid-stiffened cylinders subjected to axial compres- 
sion indicate that there is no significant difference in 
weight between optimal designs obtained with an 

Table 10 

Buckling load factors and strain level factors for grid-stiffened cylinder 
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